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ABSTRACT 



This is the second paper of a series in which we attempt to put constraints on the flattening of dark halos in disk galaxies. For this 
purpose, we observe the HI in edge-on galaxies, where it is in principle possible to measure the force field in the halo vertically 
and radially from gas layer flaring and rotation curve decomposition respectively. To calculate the force fields, we need to analyse 
the observed XV diagrams to accurately measure all three functions that describe the planar kinematics and distribution of a galaxy: 
the radial HI surface density, the rotation curve and the HI velocity dispersion. In this paper, we discuss the improvements and 
limitations of the methods previously used to measure these HI properties. We extend the constant velocity dispersion method to 
include determination of the HI velocity dispersion as a function of galactocentric radius and perform extensive tests on the quality of 
the fits. 

We will apply this 'radial decomposition XV modelling method' to our HI observations of 8 HI rich, late-type, edge-on galaxies in 
the third paper of this series. 

Key words, galaxies: structure; galaxies: kinematics and dynamics; galaxies: halos; galaxies: ISM 



1. Introduction 



In paper Brienetal I I20T0I1 in this series we presented HI 
observations of a sample of 8 edge-on, HI rich, late-type galax- 
ies. The aim of the project has been described there in detail. 
Briefly, we attempt to put constraints on the flattening of dark ha- 
los around disk galaxies by measuring the force field of the halo 
vertically from the flaring of the HI layer and radially from ro- 
tation curve decomposition. For the vertical force field we need 
to determine in these galaxies both the velocity dispersion of the 
HI gas (preferably as a function of height from the central plane 
of the disk) and the thickness of the HI layer, all of this as a func- 
tion of galactocentric radius. In addition we also need to extract 
information on the rotation curve of the galaxy and the depro- 
jected HI surface density, as a function of galactocentric radius. 
In this second paper in the series we will review earlier deter- 
minations of these properties and then describe the methods we 
have developed to analyse our sample. 

To meet our goals we need to measure the vertical structure 
of galaxies out to the lowest surface densities possible and use 
galaxies with relatively large fractions of their mass in their ha- 
los, i.e. high total mass-to-light ratios galaxies. For that purpose 
we defined a sample of nearby, HI rich, late-type edge-on galax- 
ies. From paper I we recall that for a vertically isothermal gas 
sheet with a vertically Gaussian density distribution, the gradi- 
ent of the total force in the vertical direction can be written as 



dz 



(FWHM,„/2.35)2 



(1) 
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where cr,, is the vertical velocity dispersion of the gas and 
FWHMj the gas layer thickness. 

The vertical force field comes from the gas itself, plus the 
dark halo and the stars of the disk. In regions where the halo is 
a significant contributor to K., the HI flaring will be smaller for 
the more flattened halos with a smaller axis ratio q - cja (the 
ratio of the halo polar axis c to its major axis in the galactic plane 
a). 

In the present paper we will first present a new method to 
accurately determine the rotation curve, deprojected HI surface 
density and HI velocity dispersion at all radii in a edge-on gas 
disk. The motivation for this measurement method came from 
the strong dependence in Eqn. (2) in paper I of the derived dark 
halo flattening q on the measured HI velocity dispersion and ver- 
tical HI gas disk flaring. The superposition of velocity profiles 
from many radii in each sightline through an edge-on HI disk 
tends to cause an overestimate of the velocity dispersion with 
most measurement methods, as explained below. Measuring the 
radial flaring profile requires a model of the galaxy rotation and 
face-on surface density; this also necessitates high accuracy ro- 
tation curve measurement and HI surface density deprojection. 

The HI velocity dispersion has not previously been measured 
systematically in edge-on galaxies. Most measurements of the 
HI velocity dispersion have been conducted on face-on galax- 
ies, and these measurements were largely conducted 15-20 years 
ago with relatively low spatial resolution FWHMfl= 1-3 kpc. In 
Table[T]and|2] we show the resolution of each observation used to 
measure the velocity dispersion of near face-on and near-edge- 
on galaxies, respectively. The radial structure of the HI velocity 
dispersion has been measured in only 6 galaxies - 5 from near 
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Table 1. Summary of previous observations of gas velocity dispersion in face-on galaxies. 





1 


Type 


ni veiuLiiy 


Resolution 




Hi Tcidius 




F>ict 
L.'ISL. 


Reference 










Spectral 


Spatial 





















km s"' 


kpc 


kpc 




6 


Mpc 




NGC3938 


10 ±5 


Scl 


constant*" 


8.2" 


1.2 X 1.7 


13.1 


1.7 


7.5 


10 


van der Ki'uit & Shostak (1982) 


NGC628 


5 ±5= 


Scl 


falling'' 


8.2" 


2.9 X 2.9 


40.7 


2.8 


14.0 


10 


Shostak & van der Kruit (1984) 


NGC1058 


8±2f 


ScII-III 


constant'' 


8.2" 


2.2 X 2.2"^ 


18.9 


4.3 


8.7 


10'^ 


van der Ki'uit & Shostak (1984) 


NGC1058 


8 ±2' 


ScII-III 


falling'' 


2.58" 


2.1 X2.0'' 


14.1 


3.2 


6.7 


lO"" 


Dickey et al. (1990) 


NGC6946 


34 


Scd/ScI 


falling" 


8.25" 


0.7 X 0.8 


17.2 


2.0 


21.2 


5.9'' 


Boulanger & Viallefond (1992) 


NGC5474S 


21+" 


ScdIV 


falling'' 


5.2" 


1.2 X 1.1 


10.5 


4.6 


8.6 


7 


Rownd et al. (1994) 


NGC1058 


8±2f 


ScII-III 


falling'' 


2.58" 


1.5 X 1.5 


21.8 


5.0 


15.0 


10 


Petric & RuDen (2006) 



" Hanning smoothed. 

'' Distance shown is mean distance calculated from brightest stars in NGC6946 group. 

To fac ilitate comparis on with the other ob servations for NGC1058, we have used the adopted distance of 10 Mpc used by iPetric & Rupeiil 
| |2006I) . rather than the Ivan der Kruit & Sho stak ( 1984) preferred value of 9 Mpc. 

^ To facilitate comparison with the other observations for NGC1058, we have used the adopted distance of 10 Mpc used by iPetric & Rupeiil 

(2006), rather than the Dickey et al. ( 1990) preferred value of 10.2 - 14.5 Mpc. 

Warps from a; 5° within the optical disk to i a; 0°, then back to ; s; 10° at the Hi edge ( Shostak & van der KruR ll984h . 
Warps from » 8 ± 2° within the optical disk to exactly face-on at the Hi edge (van de r Kruit & Shostak. mH) 

^ Classifi ed as peculiar, and k nown to be interacting with its much larger neighbour MlOl which is only 90 kpc distance assuming a distance of 
7 Mpc ( lRowridetallll994l) . 

'' Direct measurement of the velocity dispersion of the broadening of the vertical velocity dispersion from a near face-on disk. Measurements 
corrected for instrumental broadening and inclination, but not beam-smearing. 



face-on galaxies and one from an edge-on galaxy. The 5 face- 
on galaxies were observed at ~1 kpc or larger resolution. The 
uncertain inclination of low resolution face-on HI observations 
can cause additional broadening, due to projection of the gas 
rotation along the line-of-sight. Galaxies inclined slightly away 
from face-on will suffer beam smearing of the gradients in the 
projected rotation as well as beam smearing of the gradients of 
the intrinsic velocity dispersion. 

From Table [1] it is seen that only 3 of the low-inclination 
velocity dispersion studies were undertaken with good spec- 
tral resolution (channel FWHM,, ^ 5 km s '). Of these three 
high spectral res olution studies, only the recent measurements 
of NGC1058 by iPetric & Rupei] (l2006l) had sufficient spatial 
resolution to sample the HI radius by at least 10 synthesised 
beam FW HMs (6,„ai)- The main d ifference between the obser- 
vations of IPetric & RupenI (|2006|) and the other low inclina- 
tion HI ve locity dispersion meas urements are that the obser- 
vations by IPetric & RupenI (l2006l) were deeper (23 hr total in- 
tegration time on 3 VLA array configurations) and with good 
spectral resolution. The two low HI inclination velocity disper- 
sion studies that found a constant velocity dispersion with radius 
were only marginally spectrally resolved (channel FWHM,.= 
8.2 km s"'), and likewise also marginally spatially resolved 
(5 < Rni/FWUMg < 10). 

Of the HI velocity dispersion measurements of high inclina- 
tion g alaxies (see Table |2]i, only the study of NGC4244 (Oiling, 
1199 6) has the spatial resolution (~100 pc) to investigate the ra- 
dial velocity dispersion on scales relevant to common galactic 
structure (star formation, spiral density waves, bars). However, 
Olling's method ( 1996i) used only part of the HI emission dis- 
tribution, and may not be accurate (see Sect. 12.6b . The other 
high inclinat i on ve locity dispersion measurements obtained by 
iKregel et al.l (12004 ) used an HI velocity dispersion fit that was 
held fixed with radius, consequently no information was ob- 
tained about the radial variability of the velocity dispersion (see 
the discussion of the method in Sect. 12.3b . 

As a result, very little is known about the velocity dispersion 
of the HI in the ISM. The dependence of the HI velocity disper- 



sion on galaxy mass, surface brightness, star formation rate, and 
density variations such as spiral arms or bars is unknown. Proper 
modelling of the full HI emission cube gives the radial distri- 
bution of surface density, rotation and velocity dispersion. This 
approach could be expanded to also model the HI transparency 
of the HI emission distribution of a galaxy, whether edge-on or 
inclined. 

In the past, low telescope resolution and deprojection issues 
made such studies nonviable. In this paper, we show that it is 
possible to accurately perform the deprojection required to mea- 
sure the velocity dispersion, rotation and surface density as func- 
tions of radius in edge-on galaxies. The method is also applica- 
ble to galaxies of moderate inclination. The edge-on orientation 
of galaxies in our study is only required to minimise the uncer- 
tainty of the measured vertical gas flaring, which is required to 
determine the dark halo flattening. If fitting galaxies less inclined 
than edge-on, the inclination must be accurately fitted at all radii, 
and the galaxy must not be too close to face-on or else the line- 
of-sight velocity will be not be sufficiently resolved to model the 
galaxy rotation. Application of our new HI emission modelling 
method to nearby non-edge-on spiral galaxies, would allow in- 
vestigation of the velocity dispersion, transparency, temperature 
and density structure of neutral hydrogen in the ISM of a large 
sample of galaxies. This could explain how this small-scale ISM 
structure varies with the mass and star formation rate of spiral 
galaxies. 

Many methods have been devised to measure the rotation 
curves of edge-on galaxies. In Table |3] we display the previous 
rotation curve measurements of edge-on galaxies, showing the 
method used and the resolution of the observations. In the next 
section, we will briefly review these methods, and the few previ- 
ous methods to measure the gas velocity dispersion and surface 
density, before presenting our new method in Sect. [3] 

In the paper III of this series we will present the results of 
applying this new method to the HI observations of the edge-on 
galaxy sample presented in paper I. In that paper, we will paper 
also derive the radial flaring of the HI layer of these galaxies. A 
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Table 2. Summary of previous observations of gas velocity dispersion of highly inclined galaxies. 





1 


Type 


111 veioLiLy 


Resolution 




Hi radius 




Distance 


Ref. 










Qlspci aiUIi 


Snpptrfil 


Snfiti al 






















km s^ ^ 


kpc 


"25 


pc — 


9 


Mpc 






NGC4244 


84.5'' 


Scd 


constant 


5.2" 


0.17 


630 


II.O 


2.0 


3.6 


Oiling (1996) 


ESOI42-G24 


90'' 


Scd 


constanf^ 


6.6" 


3.8 


19.0 


1.6 


5.0 


26.6 


Kregel et al. 


(2004) 


ESOI57-GI8 


90'' 


Scd 


constant'' 


6.6 


2.6 


11. 1 


1.3 


4.3 


18.1 


Rregel et al. 


(2004) 


ESO20I-G22 


90'' 


Sc 


constant" 


6.6 


9.0 


34.1 


1.5 


3.8 


55.4 


Kregel et al. 


/ r~\f\f\ A \ 

(2004) 


ESO240-GI1 


90'' 


Sc 


constant" 


6.6 


1.7 


41.5 


1.4 


24.2 


38.1 


Kregel et al. 


(2004) 


ES0263-GI5 


90'' 


Sc 


constant" 


6.6 


5.1 


24.8 


1.5 


4.8 


35.2 


Kregel et al. 


(2004) 


ES0269-G15 


90'' 


Sc 


constant" 


u.o 


7 S 
/ .J 


19.9 


1.0 


2.7 


46.7 


Kregel et al. 


(2004) 




Qnb 
y\i 


Co 

oa 


constant" 


6.6" 


7.3 


14. J 


1 

i.Z 


z.U 




Kregel et al. 


(2004) 


ES0416-G25 


90'' 


Sb 


constant" 


6.6" 


17.7 


36.7 


1.6 


2.1 


67.0 


Kregel et al. 


(2004) 


ES0435-GI4 


90'' 


Sc 


constant" 


6.6" 


7.7 


18.9 


1.5 


2.5 


35.2 


Kregel et al. 


(2004) 


ES0435-G25 


90'' 


Sc 


constant" 


6.6" 


2.2 


36.1 


1.2 


16.3 


32.4 


Kregel et al. 


(2004) 


ESO435-G50 


90'' 


Sc 


constant" 


6.6" 


7.1 


13.1 


1.4 


1.9 


35.2 


Kregel et al. 


(2004) 


ES0446-GI8 


90'' 


Sb 


constant" 


6.66" 


12.2 


31.1 


1.3 


2.5 


64.1 


Kregel et al. 


(2004) 


ES0564-G27 


90'' 


Sc 


constant" 


6.66" 


10.2 


28.7 


1.6 


2.8 


29.5 


Kregel et al. 


(2004) 


NGC5170 


90'' 


Sc 


constant" 


6.66" 


7.4 


31.9 


1.3 


4.3 


20.9 


Kregel et al. 


(2004) 



" Hanning smoothed. 
'' From LEDA database. 

" The unwarped inner disk has an incl ination of ; = 84.5°, outside of which there is a small warp to i = 82.5 ± 1°, before warping back to the 

plane of the inner disk at large radii (Oiling, 1996). 
'' Direct measurement of the velocity dispersion of the broadening of the extreme velocity envelope of the XV map of an edge-on galaxy. 

Measurements corrected for instrumental broadening, projection and beam-smearing, via Olling's method ( see Sect. 12. 6t . 

" Iterative envelope fitting of the XV map of an high-inclination galaxy integrated over the galaxy minor axis i lKregel et al.Ll2004l) . This method 

was unable to measure radial structure of the velocity dispersion, as the method did not include beam correction and the galaxies in the sample 

were observed with very large beams. 



determination of the flattening of the dark halo of one galaxy in 
our sample, UGC7321, will be presented in paper IV. 

2. Methods to extract rotation curves for edge-on 
galaxies 

High accuracy rotation curves and deprojected HI surface densi- 
ties can be difficult to measure accurately from HI observations 
of edge-on disk galaxies, due to the projection of the velocity 
profiles from a range of radii into a single line-of-sight velocity 
profile. We will first discuss the methods that have been used pre- 
viously and the accuracy of their results. None of these methods 
independently determine the HI gas velocity dispersion from the 
data. These methods provided a useful starting point from which 
to develop a new method that incorporates the determination of 
velocity dispersion directly from the observations. 



2. 1 . The peak flux & envelope tracing method for measuring 
the rotation curve 

Most previous measurements of the HI rotation curves from 
edge-on galaxies use the shape of the outer envelope of the 
position velocity diagram, using either the Peak Flux method 
(iRubin et all 119851: iMafliewson et all 119921 lUs on & Matthews, 
2003") or the Envelope Tracing method ifSancisi & AUeni il979t 
Sofue, 1996; Garcia-Ruiz et al., 2002|). The Peak Flux method 
adopts the line-of-sight velocity at the maximum flux along each 
slice of the XV diagram as the rotation velocity at the line of 
nodes in that sightline. The Envelope Tracing method instead 
uses the line-of-sight velocity at a nominal level, e.g. 20% of the 
peak flux, on the outer velocity side of each slice through the XV 
diagram. Both these methods also involve application of a cor- 
rection to the measured velocities, to correct for the combined 



line broadening due to the instrumental broadening and the un- 
known intrinsic velocity dispersion of the HI. 

The main problem with such methods is that the flux level 
of the actual rotation velocity at any one radius is dependent 
on three different unknown radial functions: the actual rotation 
curve, the deprojected HI surface density and the HI velocity 
dispersion. The peak flux method underestimates the rotation, 
while the envelope tracing method tends to adopt the velocity at 
a flux level dependent on the signal-to-noise of the observations. 
Too low a flux level will overestimate the rotation, while the 
converse is true, if too high a level is chosen. 

The accuracy of derived rotation velocities also depends on 
the accuracy of the adopted gaseous velocity dispersion used 
in the velocity broadening correction of the measured velocity. 
Users typically adopt a nominal radially invariant value of 6 - 10 
km s ' from direct vertical HI v elocity dispersion measurements 
of normal face- on disk galaxies (Ivan der Kruit & Shostaki [19841: 
iLo et al.LfT993l) . 

Beam-smearing also tends to increase the velocity of the 
outer envelope, while lowering the velocity of the peak flux of 
slices in the inner disk regions, where the radial gradient of ro- 
tation is larger Beam smearing is a significant problem in most 
HI rotation curve studies of edge-on spiral galaxies as few mea- 
surements use observations with a resolution greater than 1-2 
kpc (see Table|3]l. 

The shape of the extreme velocity edge of the XV diagram is 
constructed from superposition of gas all along the line-of-sight, 
including gas from the line-of-nodes to the edge of the disk. 
The gas at each location has an intrinsic velocity dispersion that 
may be variable with radius. The edge-on orientation of the disk 
projects all of these radial velocity profiles onto a single line-of- 
sight velocity profile at each li ne-of-sight position through th e 
disk (see for example Fig. 1 of iKregel & van der Kruill 120041) . 
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Table 3. Summary of previous measurements of rotation curves in edge-on galaxies (/ > 80°). 



f^al a Y1 PC 


MpthnH 




Resolution 


ctfinr'P 










Spectral 


Sp3.tial 












km s ' 


kpc 


Mpc 




UGC7321 


Peak Flux" 


Hi 


5.2 


0.8 


10 


Uson & Matthews (2003) 


NGC891 


Envelope Tracing'' 


Hi 




1. / 


9.5* 


Sancisi & Allen ( 1979) 


NGC891^ 


Envelope Tracing*" 


Hi 




U. / 


n cf 

9.5 


Swaters et al. (1997) 


9 Qwaris 


Envelope tracing ■ 


Hi 


4-8 


0.5-2.8 


3.4-19 


awaters ( 1999) 


4 spirals 


Envelope Tracing'' 


HI 


on A 1 


n n o Q 


9.5-15.5 


Sofue (1996) 


4 LSBs* 


Envelope Tracing*" 


Hi 


4-8 


0.7-2.7 


5-19 


de Blok 8l Bosma (2002) 


NGC891f 


Fixed cr,, & Parametric v{R) XV 


Hi 


27 


1.7 


9.5 


van derKruit (1981) 




Modelling'''' 














Fixed (T,, & Parametric v(K) XV 


HI 


1 7 


U.o 


/.9 


DOttema et al. ( lysu) 




ivioQeuing 












NGC3079 


Iterative Envelope Tracing'' 


CO 


5.2 


0.12 


15.5 


Takamiva & Sofue (2002) 


14 spirals 


Fixed cr,, ID Gaussian 


Hi 


6.6 


1.7-17.7 


18-67 


Kregel et al. (2004) 




Envelope Fitting'' 












24 spirals 


Fixed tr,, 2D Gaussian 


Hi 


5.0-20.1 


0.8-4.8 


5-33 


Garci'a-Ruiz et al. (2002) 




Envelope Fitting'' 












8 spirals* 


Fixed cr,, XV Modelling 


Hi 


6.6-33.0 


0.9-9.8 


9.5-54* 


Kreeel & van der Kruit (2004) 


NGC4244 


Corrected Gaussian Envelope 


Hi 


5.2 


0.17 


3.6 


Oiling (1996) 



Fitting'' 



" Corrected for instrumental broadening & gas velocity dispersion using the "equivalent rectangular measure" iShane & Bieger-Smithl[T965) . 
*' Corrected for inst rumental broadening & g as velocity dispersion by quadratic subtraction, see Eqn. 

'^ Uses the standard [Fall & EfstathioiJ ( |198(1|) parametric rotation curve form, constant cr-, and a 2D exponential surface density to model the XV 
diagram. 

'* Corrected for broadening & beam-smearing (see OUing's method in Sect. 12. 6t . 
Asymmetric drift correction applied. 

* Using distance measured bv lvan derKruidfT98lhforNGC891. 

+ Re-analysis of NGC891 data observed bv lsancisi & Ailed dl979h . 

* New deeper WSRT NGC 89I data - - 144 h r. 

* Re-analysed subset of the SwatersI il999') dataset^ 

* Includes re-analysed NGC891 data observed bv iSwaters et alj ( 119971) . 



Due to this projection effect it can be inferred that the enve- 
lope tracing method would fail to accurately measure the rotation 
curve, even if it were used with a correct HI velocity dispersion, 
high-signal-to-noise and high resolution data. 

Consequently, the rotation curve derived with these two 
methods may exhibit an incoiTect shape regardless of the signal- 
to-noise and spatial resolution of the observations. For example, 
an adopted HI velocity dispersion say with an dispersion uncer- 
tainty of 3 km s"' , will cause a rotation curve uncertainty of 7 km 
s ' due to the velocity dispersion eiTor alone. As the HI velocity 
dispersion uncertainty is not known, if it is adopted from low res- 
olution face-on measurements of other spiral galaxies, the actual 
rotation curve error could be larger Consequently many high 
inclination rotation curve measurement are likely to be in error 
by between 7-20 km s"', when the effects of projection and 
beam smearing are included. Worse still, if the HI were partially 
opaque, the inner disk rotation may be further underestimated. 

2.2. The iterative enveiope tracing metliod for measuring tiie 
rotation curve 

In 2002, iTakamiva & Sofu3 (l2002l) developed the iterative en- 
velope tracing method. This improves on the standard envelope 
tracing method by automatically adjusting the initial envelope 
tracing rotation curve, until the envelope tracing fit to a model 
XV distribution built with this model rotation curve matches the 
original envelope tracing rotation curve fit to the observed XV 
distribution. To test the rotation curve fit derived with the stan- 



dard envelope tracing method, ITakamiva & Sofud (l2002h used 
this fit to build a model XV distribution which they also fitted 
with the envelope tracing method. By comparing the two fits, 
and adjusting the model rotation curve accordingly, they incre- 
mented the model rotation curve until it generated an XV dia- 
gram that produced the same envelope tracing fit as the observed 
XV diagram. The initial fit VfufiiR) to the observed XV diagram 
was used as the first model rotation curve v,„odei,i{R), with which 
the first model XV diagram XVi was built and a new enve- 
lope tracing fit Vfi, i{R) obtained. This was used to estimate a 
revised rotation curve model v„,i,dd.2{R) by comparing the new 
fit Vfit,\{R) to the original fit Vfi,fi(R). At each radius where the 
difference was greater than some tolerance Vcru, the model rota- 
tion curve was incremented by the difference Vfi,fi(R) - Vfi, \{R), 
thereby forming a revised model rotation curve Vmodei,2iR)- A 
new model XV diagram was built with the revised model ro- 
tation curve and the process was iterated, as shown below in 
Eqn. (2)-(4), until the fit v ju^niR) matched the original fit v fafi{R). 
Whereupon the model rotation curve Vmodei.ni^) was taken to be 
the best model for the observations. 

Av/,-,,„(7?) = Vfi,,0{R)-Vfi,,n(R) (2) 

(R) = VfiM + Avfi,,„(R) (3) 

if Avfi,^„(R) > 

(R) = v/,v,o(/?) (4) 

Takamiya & Sofue employed several assumptions to build 
each model XV diagram. Firstly, the radial gas surface den- 
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Fig. 1. Diagram of galaxy coordinate system used to model the 
XV gas distribution in Eqn. (|6]). This shows the major axis /?', 
depth along the line-of-sight s, and the projection of galaxy ro- 
tation along the line-of-sight towards the observer 



sity was approximated by the gas surface density along the ma- 
jor axis, which was derived from the observed XV diagram 
by integrating over velocity. Secondly, the gas velocity dis- 
persion was assumed to be isotropic and isothermal in R and 
z with an adopted a value of 10 km s which is typically 
observed in HI in spiral disks (Shostak & van der Kruit, 1984; 
Ivan der Kruit & ShostakLlTosl . In the z direction the HI surface 
density was approximated by a sech^ model with a constant scale 
height of 100 pc. To simulate the observing process, each model 
XV cube was then convolved by the spatial and spectral resolu- 
tion of the telescope, and integrated over z to form the model XV 
map. 

iTakamiva & Sofu3 (12002 ') applied this method and the ear- 
lier peak flux and envelope tracing methods to optical emis- 
sion spectra and CO (J=l-0) observations with the Nobeyama 
Millimetre array on the edge-on galaxy NGC3079, and four 
moderately inclined galaxies: NGC4536, NGC4254, NGC4419, 
NGC4501. They found that both the envelope tracing and peak 
flux methods underestimated the rotation velocities in the in- 
ner regions of their galaxies, with the peak flux method be- 
ing in error by up to 50 - 100%. Furthermore at outer radii 
with low signal-to-noise, the derived rotation velocities ex- 
hibited dramatic scatter of +20 - 100 km s"' (see Fig. 4 of 
[T akami ya & Sofue, 2002), probably because the 3cr noise level 
was above the 20% flux level used for envelope tracing. 

[Sofue et al.l (l2003h applied the method to CO observations 
of a sample of moderately inclined galaxies (25° < ; < 74°), 
but this time with a small difference. Instead of fitting the XV 
diagrams at a single fractional flux level, they fitted the line-of- 
sight velocity at a series of flux levels /= 20, 40, 60, 80 and 
100%, and replaced Eqn. (|2]l by 

Av/,v,„(/?) = X ^' ^ [4,o(^) - (5) 



where gi is a set of weights satisfying YjI 8i- 1> ^nd usually 
evenly weighted, except when low flux levels are biased by 
noise in sightlines through the XV diagram with low signal-to- 
noise. The iterative envelope tracing method is an improvement 
because it checks the veracity of the measured rotation curve. 
However, substantial flaws remain: it uses an incorrect radial gas 
surface density, and an assumed constant gas velocity dispersion. 

2.3. The fixed velocity dispersion, 1 D Gaussian envelope 
fitting method for measuring the rotation curve & 
approximating the gas velocity dispersion 

iGarcfa-Ruiz et al.l (12002 ') improved on the standard Gaussian 
model of the full outer envelope of the XV diagram, by perform- 
ing a series of three-parameter Gaussian fits to the outer envelope 
of the velocity profile in each sightline. The first channel on the 
outer velocity envelope with flux > 3cr was the outermost ve- 
locity channel used in each Gaussian fit. For the first Gaussian 
fit, this outer channel and the two lower channels were included 
in the profile sample. Subsequent fits were performed, each time 
adding the next lower velocity channel, until the velocity chan- 
nel with peak flux was reached. The centre of the Gaussian with 
the smallest reduced was retained for that sightline. 

iKregel et al.l (|2004|) applied the fixed velocity dispersion ID 
Gaussian envelope fitting method to 14 galaxies using an radi- 
ally constant velocity dispersion of 10 km s Successful fits 
were obtained for 11 galaxies using this HI velocity disper- 
sion, however the remaining three galaxies required the veloc- 
ity dispersion to be reduced to 6 - 7 km s"' to achieve a con- 
verging fit. However, the large beam size of the observations 
(1.7 i FWHMfl i 18 kpc) masked any small scale variations in 
the rotation speed. Only 3 of the galaxies in their sample have 
a HI radius greater than 5 beamwidths. (see Table |2]i. For two 
of the three galaxies, which required smaller HI velocity disper- 
sions, the radius was spanned by only 2.5 synthesized beams. 
Consequently, due to the low spatial resolution, they were un- 
able to measure any radial structure of the velocity dispersion. 

We show in Sect. I4.3l that even an intrinsically constant gas 
velocity dispersion should appear broadened near the Galaxy 
centre due to beam smearing — typically by i 50 - 100%, if 
FWHMg > 1 kpc. This beam smearing eff'ect would cause the;^^^ 
error of the velocity dispersion fits at inner radii to dramatically 
increase as the XV slice bears progressively less resemblance to 
a sum of Gaussians. 

An advantage of this method is that it could be extended to 
measure the gas velocity dispersion as a function of radius, by 
allowing the velocity dispersion to be a free parameter in the 
Gaussian envelope fits, and correcting the measured velocity dis- 
persion for instrumental broadening using Eqn. (|9]l. 

2.4. The fixed velocity dispersion, 2D Gaussian envelope 
fitting method for measuring the rotation curve 

To attempt to correct the fixed vel ocity dispersion ID Gaus sian 
envelope fitting method (Sect. l2.3l l. lGarcia-Ruiz et aT1(l2002h ap- 
plied the same series of fixed dispersion Gaussian fits to slices 
along the radial R' direction of the XV diagram; i.e., in addi- 
tion to fitting at velocity profile at each major axis position R', 
they also fitted major axis R' profiles of the flux in each veloc- 
ity channel v. The width of the outer envelope in the R' direc- 
tion was determined by the FWHMe of the synthesised beam 
and the unknown gradient of the rotation curve. They measured 
this spatial dispersion along the R' profile at the systemic ve- 
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Fig. 2. This is a derived gas kinematics fitted from a raw synthetic XV map (ie. not spatially smeared by the telescope beam and 
without noise) using the radial decomposition method. The plots of gas velocity dispersion, rotation and deprojected surface density 
show the actual synthetic galaxy kinematic as a grey (yellow) line and measured fits in (red) crosses. The plot in the bottom right 
hand corner shows both sides of the measured rotation curve (in the color version the receding side is plotted first in red and the 
approaching side on top in blue). Error bars are shown on all fits, though in most cases they are smaller than the plot symbol. 
Subsequent similar figures in this paper will use the same plotting scheme. 



locity, vios=Vsys, of each galaxy, and held this value fixed when 
Gaussian fitting the R' profile of each channel. The rotation ve- 
locity at each radius was taken to be the minimum of the rotation 
velocities measured with each method. . 

Usi ng an adopted HI velocity dispersion of 8 km s 
iGarcfa-Ruiz et alJ (l2002h applied this method to measure the ro- 
tation curves of 24 edge-on (/ > 80°) galaxies (see Table [3]). 
Beam-smearing correction was important for this sample as most 
of their observations had a FWHM beam of 1 -2 kpc, and most of 
the galaxies in their sample were moderately resolved spatially 
with 5 < Rni/FWUMg < 10. 



tion curve and the radial HI distribution. It was first ap- 
plied by van der Kruit (1981) and was further developed by 
iKregel & van der Kruit ( 20 0J). We regard this method as the 
most powerful approach. Kregel & van der Krui i (l200l model 
the full observed XV diagram as a set of concentric, coplanar 
rings, by progressively measuring the rotation v(R) from the out- 
side in, or as they call it the "onion-peeling" approach. 

For a diagram of the coordinate system we refer to Fig.[T] By 
modelling the HI layer as an axisymmetric, transparent isother- 
mal sheet in pure circular rotation, the velocity profile probed by 
sight lines through the XV diagram is 



2.5. The fixed velocity dispersion, XV modelling method for 

measuring the rotation curve and deprojected HI surface r s _ 1 f ds I. (R) (6) 
density ' ct,,^, ^ J 

The basic approach of this method is to directly fit the full / [y - vdR) (R' /R)]^ 

observed distribution in the XV plane to derive the rota- ^^Py "^crlto, 
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Fig. 3. Flow chart showing the iterative radial decomposition 
method. 



where 

V/, 



(7) 



and R' is the projected radius on the sky, s the depth along the 
line of sight (i.e. s- at the line-of-nodes), R the galactocen- 
tric radius, 'Lhi(R) the radial face-on surface density profile and 
o'v,tot(R) the total velocity dispersion; see Eqn. (|9]i. vdR) is the 
rotation velocity and vios{R) is the projection of the rotation ve- 
locity onto the line-of-sight. This method assumes a constant 
HI velocity dispersion cr^.jfj of 10 km s"'. The face-on surface 
density UniiR) was derived independently using the Warmels' 
119881) meth od, which uses the Lucy alg orithm to invert an Abel 
integral (e.g.'Bin nev & Tremainelll98 7. Eqn. (IB) - (59b)). 

At each sightline the HI surface density and rotation on the 
line of nodes was fitted using a single Gaussian. By starting at 
the outermost projected radius 7?' =7?^,,,^^ and using at least 4 in- 
dependent data points above a nominal noise flux level, the ro- 
tation velocity was determined by a Gaussian fit to the outer en- 
velope. The next inwards profile at /?'=(/?^„,p^ - AR') was then 
fitted simultaneously for {Router ~ ™d R' ^ where tsR' de- 
noted the annulus width, which was set to the pixel size of the 
observations. An initial estimate of the value of the rotational ve- 
locity determined in the previous fit was used at R'ouier^ while at 
the line-of-nodes (/?^„,^,.- A/?') the ID Gaussian envelope tracing 
value was measured. Each line profile was fitted as a set of fixed 
dispersion Gaussians with flux determined from the derived sur- 



face density, yielding the rotation velocity at the line-of-nodes 
v{R - R') and new values for the rotation velocities (measured 
in projection) of all high radii annuli in the sight line. 

The procedure was repeated, progressively determining the 
rotation velocity at all radii inwards until the centre was reached, 
and then the same method was applied to the other side of the 
galaxy. In this way the fitting algorithm does not become ill- 
constrained, as the number of free parameters increases at the 
same rate as the initial estimates of the fit parameters. 

In addition to more accurately determining the rotation ve- 
locity by modellin g the full profile through e ach sight line of 
the XV diagram, Kregel & van der Kruit (2004) noted that this 
method also successfully corrects the rotation velocity for the ef- 
fects of beam smearing. This was done by first building a model 
mini XV diagram comprising sufficient sightlines adjacent and 
outwards of R' to span at least the beam FWHMg using Eqn. (|6]l, 
and convolving by the instrumental broadening and any Hanning 
taper (if applied). The mini XV diagram is then convolved spa- 
tially by the known beam FWHMg. The velocity profile of in- 
terest, can then be corrected for the modelled beam smearing 
effects determined from the mini XV diagram, prior to Gaussian 
fitting. 

The best fit rotation velocity at each radius w as then de- 
termi ned by minimization using the 'amoeba' (Press et all 
1992 ) implementation of the simplex method ( Nelder & Meag, 
1965 '). This easily handles a large number of free parameters 
and does not require the derivatives of the model with respect to 
those parameters (a computation-intensive calculation). Errors 
were estimated using Monte Carlo simulations. 



'Kregel & van der Krui3 ([2004) applied this method to the HI 
observations o f 14 edge-o n galaxies. These comprised 13 of the 
galaxies in the iKregel et al. (2004) sample and HI observations 
of NGC891 bv ISwaters et al., (,1997.) . The fitting method suc- 
cessfully converged to sensible results for all the galaxies, ex- 
cept five galaxies with low peak signal-to-noise of S IN < 15, 
and one galaxy NGC5529 which displays a pronounced warp 
of the disk plane. The resolution of the galaxies with success- 
ful fits are shown in Table [3] Despite the large beam sizes 
(1 i FWHMfl i 10), two of these successfully fitted galaxies 
(NGC891 and ES0435-G25) that were well sampled by the 
beam with /^////FWHMg > 10, while four had moderate sam- 
pling with 5 < 7?////FWHMe < 10, and two were successfully 
fitted despite poor sampling of /?////FWHMe < 5. 

For four galaxie s with higher S/N observations, 
iKregel & van der Kruitj (l2004l) also succeeded in directly 
fitting the surface density by allowing it to be a free parameter in 
the Gaussian fits. They compared the rotation curve and surface 
densities derived from their method with the rotation curve 
derived from standard envelope tracing method and the surface 
density derived from the Warmels' (1988) deprojection method. 
Each method was applied to a simulated HI XV diagram of 
NGC2403, built with the "true" rotation curve and surface 
density measured by Fraternali et al. (2002) using extensive 
3-D modelling, showing that the fixed velocity dispersion 
XV modelling method correctly measures both the rotation 
curve and the deprojected surface density. Kregel & van der 
Kruit's surface density deprojection method is a significant 
improvement on that of (Warmels, 1988), which was unable to 
recover the surface density variations on radial scales smaller 
than the FWHMg of the synthesised beam. 
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Fig. 4a. Actual and observed XV maps of the synthetic galaxy. This has the properties illustrated by the grey (yellow) lines in 
Figs. |4b] and |4c] The top XV diagrams show the actual XV HI distribution, while the lower plot is the XV diagram as observed 
to detection limits corresponding to peak noise-to-noise {S IN) of 80 and 20, each with the same telescope FWHM of 300 pc. The 
contours in the upper panels are (in units of 10'^ atoms cm ' pixel ') at 0.5, 1, 3, 10, 30, 100, 300 and 1000 and in the lower panels 
at 2, 3, 5, 7, 10, 15, 20, 30, 50 and 70 times the rm.s. noise in the channel maps. The effect of increased noise shows clearly the 
reduced fraction of the actual XV diagram that is detected. The vertical dashed lines on the observed XV maps mark the outermost 
radius with sufficient flux above 2cr to allow kinematic fitting. Note that the scales are different in the plots on the right and the left, 
because of changing observable extent for different S/N. 



2.6. Direct outer envelope fit metliod for measuring tine 
veiocity dispersion, rotation curve and deprojected HI 
surface density 



lOUind d 19961) used the simple outer envelope velocity disper- 
sion fit, but then endeavoured to correct the measured broaden- 
ing for the major effects of beam-smearing and velocity profile 
projection in addition to instrument al broad ening. His correction 
method (detailed in Appendix A of 'Oiling', 7996) calculated the 
effect of the projection given in Eqn. (|6]l to determine how to 
correct the measured direct Gaussian fit to the outer velocity en- 
velope of each slice of the XV diagram. By approximating the 
Gaussian in Eqn. (|6l) as a parabola, and integrating over several 
velocity channels around the line-of-node {R=R'), Oiling derived 
an analytical approximation of t he velocity pr ofile around the 
line-of-nodes to be (Eqn. (A5) of I011ind.ll996h : 



I(R',v) 



g 2\£r„(»); X 

o-,(R) V2^ 

elW(R)v,(R) R'ei 



2- 



3cT,{R)^ 



aW(R) , / W(Rf 
+ b 



c 



cr.iRV 



1 + 



(8) 



where the kinematics (ie. v^R), cr,,(/?) and 'Eh[{R)) were evalu- 
ated at the line-of-nodes R-R', and W(R)-v - vdR). Av is the 
channel width, K, a, b and c are constants, em-SmaxIR' , and s is 
the depth along the line-of-sight (see the schematic of the coor- 
dinate system in Fig. [T]l. The first term, 2, results from gas on 
the line-of-nodes. The second term arises due to material close 
to the major axis, and the last 3 terms arise due to radial gradi- 



ents in the rotation velocity, gas velocity dispersion and surface 
density, respectively. 

[piling (1996) applied this method to one galaxy, NGC4244, 
which was observed at high spatial and spectral resolution with 
a FWHMfl of 170 pc and a FWHM channel width of 5.2 km s ' . 
Although NGC4244 is a small Scd galaxy with maximum HI ra- 
dius of 1 1 kpc, the high spatial resolution allowed the radial HI 
structure to be sampled by 63 FWHM beams on either side of the 
galaxy centre (see Table|2]i. By testing the above method on XV 
simulations of NGC4244, he found that Eqn. (jSj was sufficient 
to recover the kinematics in the outer galaxy. However, within 6 
kpc from the galaxy centre, beam-smearing became a problem. 
To derive the ga s properties in the inner disk consistent with the 
beam smearing. lOUingl d 19961) used simulations of NGC4244 to 
determine an empirical correction to the kinematics measured 
with Eqn. (jSj. In order to simulate the XV diagram of NGC4244, 
Oiling had to estimate the intrinsic kinematics, as the derived 
corrections depend on the intrinsic XV diagram as well as the 
resolution of the observations. With this method Oiling found 
that the velocity dispersion of NGC4244 was roughly constant 
at 8.5 + 2 km s"'. Unlike the HI velocity dispersion measure- 
ments of face-on galaxies which mostly decline with radius (see 
Table [U, Olling's analysis of NGC4244 shows no radial struc- 
ture in the velocity dispersion despite being measured on scales 
of -100 pc. 



3. The radial decomposition XV modelling method 

We developed the radial decomposition XV modelling method 
to more accurately measure the rotation, velocity disper- 
sion and deprojected surface density of HI disks. Like the 
Kregel & van der Kruit (2004) method it operates by fitting the 
galaxy from the outer radii inwards - "onion-peeling" as they 
called it. The radial decomposition method works by progres- 
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Fig. 4b. Fitted kinematics measured from an XV diagram observed with a telescope beam of FWHMe=300 pc and to a peak signal- 
to-noise of 80. Beam-smearing correction is not necessary at this high resolution, as the scatter due to beam-smearing is much 
smaller than the scatter due to noise. This figure and the plots in Fig.|4c]demonstrate the effect of noise on the measured kinematics. 
To reduce the scatter we have averaged the measurements using a evenly weighted 5 pixel wide bin after fitting the whole XV 
diagram via the radial decomposition method. The outermost dashed lines indicate the maximum radius R„,ax with sufficient flux to 
be fitted. The other dashed vertical lines mark a beam HWHM from R,„ax on either side, and a HWHM from the galaxy centre. 



sively fitting each sightline of the observed XV diagram, starting 
at the outermost edge, to measure the rotation velocity, surface 
density and ve locity dispersion. The basic me thod is essentially 
the same as in lKregel & van der KruitI (|2004|) . however we also 
fit the HI velocity dispersion, and extend the method to iterate 
further if required. By progressing inwards, the method uses the 
measured kinematics of outer radii to isolate and fit only the ve- 
locity profile of the line of nodes gas along each sightline. 

This allows one to accurately fit the rotation, velocity dis- 
persion and surface density as functions of galactocentric radius 
despite the edge-on orientation of the galaxy. The measured ve- 
locity dispersion is corrected for the broadening due to the finite 
spectral resolution and the derived HI velocity dispersion at each 
radius is used in the measurement of all inner annuli. 

The method works by modelling the radial and velocity dis- 
tribution of each slice through an edge-on galaxy that has been 
integrated over the vertical axis z- This requires the usual as- 
sumptions of circular motion, azimuthal symmetry, isotropy of 



the HI velocity dispersion tensor. To maximise the peak signal- 
to-noise of the XV map, the map is formed by integrating over 
the full range z spanned by the edge-on galaxy in the vertical 
direction. This requires that the velocity dispersion is vertically 
isothermal and that the rotation is cylindrical. In future studies, 
the assumption of HI isothermality should be tested by applying 
the radial decomposition XV modelling method to deep observa- 
tions of both low and high lattitude HI observed at high spatial 
resolution. This would determine whether the velocity disper- 
sion (and rotation) varies with height above the galactic plane. 
Althoug h the HI in some deeply observed galaxies like NGC 
891 (■Fraternahetalil2005h does show a rotational lag in the HI 
layer a few kpc above the galactic plane, our data is for the most 
part confined to a thin layer of FWHM < 2 kpc, so the weighted 
effect of the high-latitude lag on our analysis is unlikely to be 
significant. 

Although the method will not detect spatial asymmetries in 
the gas surface density, such as might be expected from bars or 
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Fig. 4c. Fitted kinematics measured from an XV diagram observed with a telescope beam of FWHMe=300 pc and to a peak signal- 
to-noise of 20. Beam-smearing correction is not necessary at this high resolution, as the scatter due to beam-smearing is much 
smaller than the scatter due to noise. This figure and the plots in Fig.|4b]demonstrate the effect of noise on the measured kinematics. 
To reduce the scatter we have averaged the measurements using a evenly weighted 7 pixel wide bin after fitting the whole XV 
diagram via the radial decomposition method. The dashed line are as in Fig.l4bl 



spiral waves, it will detect radial annuU devoid of HI, or a central 
Hi hole, as the surface density is modelled as 

The method fits the radial kinematics by modelling each 
sightline through the XV diagram as a cut through a superpo- 
sition of annular rings, in which the gas at each radius exhibits a 
Gaussian velocity profile described by the radial gas kinematics 
seen in projection according to Eqn. Each slice comprises 
gas at a range of radial bins, with the gas in each radial bin in the 
slice forming a Gaussian profile 

• centred at the rotation velocity seen in projection, v/oj 
{vios-v{R) along the line-of-nodes, vio., < v(R) for radii outside 
the nodal radius); 

• broadened to a velocity FWHM due to the effect of intrinsic 
velocity dispersion and the instrumental resolution, FWHM,yor 
and 

• with the integrated flux of each Gaussian consisting of the 
flux at that radius. 

Starting at the outermost radius on one side of the galac- 
tic centre, each line-of-sight is modelled as a sum of Gaussians, 



where each Gaussian simulates the gas at that radius. Due to the 
approximately exponential outer Hi distribution of most galax- 
ies, the gas from radii near the line of nodes dominates the flux 
in each velocity profile slice through the XV diagram. At the 
outermost line-of-sight, the velocity profile is approximated as a 
single Gaussian as the flux can be approximated as coming from 
a single radial bin with R'-Rmax- The properties v{R), <jy Hi{R), 
l.HiiR) at that radius can be directly measured at 7?,„f„. In the 
next innermost radius, R2=Rmax - (Ipixel), the observed velocity 
profile is modelled as a superposition of the Gaussian with the 
known R^ax kinematics, and the Gaussian with kinematics from 
the line-of-nodes at radius which is fitted for rotation, velocity 
dispersion and total flux. By proceeding inwards in this manner, 
the model never becomes ill-constrained as there are only 3 free 
parameters in the fit to any sightline. 

The radial kinematics are thus obtained with radial sampling 
corresponding to the pixel size of the Hi image cube along the 
galaxy major axis R'. The derived kinematics from fitting to all 
slices of the raw XV diagram using in the above slice fit are 



J.C. O'Brien et al.: The dark matter halo shape of edge-on disk galaxies II 



11 





Fig. 5a. Fitted kinematics (without beam correction) measured from a noiseless XV diagram observed with a telescope beam of 
FWHMe=l kpc. The plots shows that the rotation velocity exhibits a small beam-smearing effect and the velocity dispersion is 
broadened due to beam-smearing by ~5Q% at the Galaxy centre. The dashed line are as in Fig.l4bl 



shown in Fig.|2] The fit has been made to a synthetic galaxy with 
a constant velocity dispersion, a simple rotation curve and an ex- 
ponential surface density distribution. The plot demonstrates that 
the radial decomposition method does an excellent job of recov- 
ering the actual galaxy structure and kinematics. Near the edge 
of the edge-on galaxy, the fits are systematically high in both the 
velocity dispersion and the surface density, because the model 
assumes that the outermost detected slice contains flux at a single 
radial bin. The overestimation of the velocity dispersion, causes 
the fitting process to compensate on the next few inner radii slice 
fits; thus the velocity dispersion fit oscillates around the actual 
value in the outer few slices. The small deviation in the centre 
slice is a numerical artifact. We remove this artifact by exclud- 
ing the central fits within 0.5 kpc of the galaxy centre. Although 
these plots show the derived fits to the XV diagram formed from 
a simple galaxy model with a constant velocity dispersion, sim- 
ple surface density and rotation curve, the same fitting success 
is achieved for fits to synthetic galaxies with more complicated 
intrinsic kinematics. 



3.1. Instrumental broadening correction 

Fits of real galaxy observations also require correction for the 
instrumental broadening due to the signal response of the spec- 
tral channels of the interferometer correlator For simplicity, 
the signal response function of the correlator channels is ap- 
proximated as a Gaussian. Consequently, the dispersion due to 
both the HI velocity dispersion and the instrumental broaden- 
ing is the quadratic sum of each property. Thus, the measured 
velocity dispersion is related to the intrinsic HI velocity dis- 
persion and the dispersion cr,, ,,,,,;. of a spectral channel with 
FWHM,, ,„,,r=2cr^, ,„,,r V21og2 by: 

^^:"'< = V'^^.m + (9) 

This allows the intrinsic HI velocity dispersion to be easily de- 
rived from each Gaussian fit. 

3.2. Beam correction 

Any method that aims to recover spatial resolution on scales of 
the image pixels must also correct for the convolution effect of 
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Fig. 5b. Fitted kinematics (without beam correction) measured from a noiseless XV diagram observed with a telescope beam of 
FWHMe=2.0 kpc. The plots shows that beam-smearing by a 2 kpc beam causes a large error in the inner rotation velocity and a 
~70% increase of the measured central velocity dispersion. The dashed line are as in Fig.l4bl 



the telescope synthesised beam. The beam is approximated an 
elliptical Gaussian, which is near circular for objects at an appro- 
priate declination relative to the latitude of the observatory. The 
telescope beam causes a smearing effect along the spatial axes 
of the Hi channel map cube by convolving the flux distribution 
with a 2D Gaussian. As the XV map is formed by integrating 
over z, then providing the beam is near circular, the beam effec- 
tively smears the spatial axis of the XV map by a ID Gaussian 
with FWHM equal to that of the beam, along the R' axis. 

Beam smearing is a problem when the galaxy kinematics 
vary strongly on intervals along the major axis smaller than the 
FWHM of the synthesised beam. In practise, the largest varia- 
tions usually occur in the inner region where the rotation curve 
rises steeply and to a lesser extent in any region with variable 
kinematics due to heating processes in the ISM. 

In the following section ( 13.31 ), we describe in detail how we 
incorporate an iterative process to check the derived kinemat- 
ics measured from the observed XV map. Briefly the process 
works as follows. Firstly, the first-pass kinematics from the the 
kinematic model is used to build a synthetic galaxy XV diagram 
using Eqn. (|6]l, which is then convolved by the known telescope 



resolution to simulate the actual observational effects. This sim- 
ulated XV diagram is then fitted with the radial decomposition 
method and the resulting kinematics compared to the original fit- 
ted kinematics. Where the difference between the kinematic fits 
is greater than a critical level the kinematic models are incre- 
mented and the process repeated until the fitted kinematics and 
the simulated galaxy resemble the original fitted kinematics and 
the observed galaxy distribution. 

Initially beam correction was attempted by iterating the fit- 
ting process as described above. However the velocity disper- 
sion fits are particularly sensitive to noise, and it was found to 
be difficult to correct the velocity dispersion by iterative fitting 
if the first pass fits were too deviant from the true fit, as of- 
ten occurs when fitting galaxies with steep inner rotation curves 
that have been observed with a large beam. The optimal method 
was found to be using the be am correction method applied by 
iKregel & van der KruitI (l2004ft whereby instead of fitting the ac- 
tual observed slice, the slice is first corrected by subtracting 
the flux contamination due to beam-smearing of flux in nearby 
slices. The residual profile is then fitted using the radial decom- 
position method. As the radial decomposition method fits pro- 
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Fig. 6a. Fitted kinematics (with beam correction) measured from XV diagram observed with a telescope beam of FWHMfl=l kpc 
and with a peak signal-to-noise of 80. To reduce the scatter in the measured kinematics, the velocity dispersion and rotation curve 
were averaged (post fitting) over an evenly weighted bin of width of 5 pixels. In Fig. |6b] the XV diagram built with the best-fit 
measured kinematics and observed to the same telescope specifications is compared to the original XV diagram. 



gressively inwards, the kinematics are known at radii larger than 
the slice of interest R', thus it is possible to form a raw "mini-XV 
diagram" comprising solely the slices at higher R'. However to 
correctly derive the flux contamination due to beam smearing a 
raw "mini-XV diagram" must be built spanning sufficient slices 
on both sides, but excluding the slice of interest. This mini-XV 
diagram is convolved along R' by a Gaussian of FWHM equal 
to that of the beam FWHMg along R', to determine the beam- 
smearing component of the flux in the slice of interest. 

Consider the observed beam-smeared velocity profile at the 
slice of interest Fq. This velocity profile is the sum of the in- 
trinsic profile /o weighted by the beam centre bo and the flux 
of all nearby slices weighted by their respective beam weights 
b„, where b„ becomes negligible for slices further than a few 
FWHMfl's from the slice of interest. Thus 



^'O = fobo + ^ fnb„ + ^ fnb„. 



(10) 



Here each beam weight b,, is the flux of a normalized Gaussian 
in the n''' pixel from the beam centre. Therefore, the intrinsic 
velocity profile is 



^0 - (Zr=l fnbn + Z„Zl fnbn) 

bo 



where bo is the peak flux of a normalised Gaussian beam 



bo^ 



1 



V27rcrg 



with (Tg - 



FWHM, 







(11) 



(12) 



Beam correction is thus undertaken by determining 
^'n=\ fnbn and 2n^i fnbn by convolving a raw mini-XV diagram 
formed from all slices but the slice of interest. Then this beam 
smearing component is subtracted from the observed slice, and 
the normalisation bo is applied. The raw un-beam-smeared mini- 
XV diagram is built using Eqn. This requires the intrinsic 
kinematics which are only known for radii atR > R' determined 
from fits outside the slice of interest. To estimate both parts of the 
mini-XV diagram, we conduct an initial radial decomposition fit 
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Fig. 6b. Top: Original XV diagram observed with a FWHMfl=l kpc telescope beam and peak 5 /A^ of 80. Middle: The XV diagram 
built using the best-fit measured kinematics and observed with the same telescope resolution and sensitivity limits. The measured 
kinematics used to build the middle plot are shown in Fig.|6a] Bottom: The difference between the top and middle diagrams. Contour 
units are 2, 3, 5, 7, 10, 15, 20, 30, 50, 70 times the r.m.s. noise in the channel maps. 



across the whole observed XV diagram. This first loop fit is ac- 
curate, except near the galaxy centre, where the high gradients of 
rotation and velocity dispersion over R' distort the line profiles. 
The kinematics near the galaxy centre are estimated (for beam 
smearing correction) by extrapolating inwards using the higher 
quality fits over the rest of the galaxy. Subsequent iterations of 
the radial decomposition fitting process can then be beam cor- 
rected using the mini-XV diagram method, whereby the mini- 
XV diagram is formed using these estimated kinematics. 

This beam correction is essentially a deconvolution of the 
XV map to remove the effect of the telescope beam. As such the 
sampling of the beam, i.e. the number of pixels per beam, affects 
the success of the deconvolution, as the deconvolution process 
amplifies the noise. The recommended sampling for deconvolu- 
tion is 3 pixels per beam FWHMe. To recover a slightly higher 
spatial resolution, which is useful when modelling the more dis- 
tant galaxies, we used a beam sampling of 5 pixels when forming 
the Hi image cubes. Both sampling densities work well, how- 
ever a yet higher sampling density ( i 8) dramatically reduces 
the peak signal-to-noise of the beam corrected XV slice, as the 
telescope beam naturally limits the spatial resolution of the ob- 
servations. As the galaxies in our sample lie at recession veloc- 
ities from 400 - 2500 km s"', and distances of 3 - 30 Mpc, the 
FWHMe beam size of the galaxy observations is 0.1 - 1.0 kpc. 
Given that the kinematics can be recovered on a scale of 1/5 of 
the beam FWHMg, this leads to a recovered spatial resolution of 
20 - 200 pc for our galaxies. 



The beam correction is necessary for galaxies with large 
beam sizes t 500 pc, as for these galaxies the beam smear- 
ing is sufficient to distort the slice shape away from a sum of 
Gaussians, over the major axis range where the kinematics vary 
significantly over the span of the beam. Without beam correc- 
tion, radial decomposition fitting of the inner disk will tend to 
oscillate with a large scatter, due to the intrinsicly non-multi- 
Gaussian shape of the velocity profiles in these sightlines. This 
occurs as an erroneous fit to gas on the line of nodes in any one 
sightline will cause the fit to the next inwards radius to overcom- 
pensate, resulting in an oscillating fit over the beam-smeared do- 
main. 

3.3. Iterating to ensure self-consistent fits 

To assess and improve the quality of the fitted kinemat- 
ics, the process is iterated by building a model galaxy 
from the original fitted kinematics, simulating the obser- 
vations and comparing the fitted properties of the model 
fnlv{R), cTy HiiR), ^Hi{R)] galaxy XV diagram to the original fit- 
ted kinematics /i [v(/?), o",, A schematic showing 
the method is shown in Fig. [3] In the initial loop, loopo, the ob- 
served XV map is fitted to obtain a rough kinematic model that is 
used to perform beam correction. This mini-XV diagram is then 
used in the next fit to the observed XV map, which we will call 
loop\. The radial functions of Hi surface density, rotation and 
velocity dispersion fitted in loop\ fit, which are the first fits in- 



J.C. O'Brien et al.: The dark matter halo shape of edge-on disk galaxies II 



15 





Fig. 7a. Fitted kinematics (with beam correction) measured from XV diagram observed with a telescope beam of FWHMfl= 1 kpc 
and with a peak signal-to-noise of 30. Although the rotation velocity and the surface density are well recovered, the plot of the 
measured HI velocity dispersion shows that a peak signal-to-noise of 30 is not high enough to perform beam-smearing correction 
for this kinematic model and moderate spatial resolution. The measurements are binned over 3 pixels post-fitting, however averaging 
only improves the measurement when it is accurate. By comparison, the observation with a larger beam to the same signal-to-noise 
recovers more of the HI emission distribution and thus recovers the radial shape of the HI velocity dispersion (compare to Fig.|9a]i. In 
Fig. |7b] the XV diagram built with the best-fit measured kinematics and observed to the same telescope specifications is compared 
to the original XV diagram. 



corporating beam correction, are then used as the first kinematic 
model for a new XV diagram. By iterating the fitting process, 
and incrementing the kinematic model between fits, the process 
is repeated until the fitted kinematics from the model equal the 
original fitted kinematics with beam correction to within a criti- 
cal threshold. 



By intrinsically incorporating fitting of the gas velocity dis- 
persion into the radial decomposition method, it is possible to 
replicate the observed XV diagram from the best fit kinematics 
as all three functions characterising the XV map are obtained: 
rotation curve v(R), HI velocity dispersion o-y [fj(R) and the de- 
projected HI surface density S///(/?). Unlike earlier methods of 
measuring rotation curves or HI velocity dispersion, this method 
is internally self-consistent. 



4. Quality testing on synthetic galaxy observations 

The quality of the derived kinematics is quantified in two ways. 
Firstly, we calculate the total absolute difference between the 
best fit and original fit for each of the three functions: v(R), 
Cv,Hii^) ^iid ^Hi{R)- This difference is summed over all mea- 
sured slices across the galaxy and divided by the number of mea- 
sured slices to present a average difference at each radius. 

Secondly, the XV map formed from the best fit kinematics 
is compared to the observed XV map by computing the resid- 
ual. The residual map (with absolute values for the residuals) is 
normalised by the observed XV map to obtain a total percentage 
difference between the actual galaxy XV diagram and the best 
fit model. The residual XV map is a good indicator of surface 
density and rotation curve accuracy, but is relatively insensitive 
to errors in the measured velocity dispersion. 
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Fig. 7b. Top: Original XV diagram observed with a FWHMg^ 1 kpc telescope beam and peak S /N of 30. Middle: The XV diagram 
built using the best-fit measured kinematics and observed with the same telescope resolution and sensitivity limits. The measured 
kinematics used to build the middle plot are shown in Fig.|7a] Bottom: The difference between the top and middle diagrams. Contour 
units are 2, 3, 5, 7, 10, 15, 20, 30 times the rm.s. noise in the channel maps. 



4. 1 . Effect of noise 

The noise level is of most concern when fitting observations that 
suffer significant beam-smearing, as the beam correction method 
uses deconvolution, which is particularly sensitive to noise. To 
isolate the effect of noise from that due to beam-smearing, we 
tested the radial decomposition method on high spatial resolu- 
tion XV maps (observed with FWHMe= 300 pc) with different 
sensitivities corresponding to peak signal-to-noise ratios rang- 
ing from 20 to 80. In each case the same galaxy model is used, 
featuring a differential rotation curve, an exponential face-on HI 
surface density (with a HI scalelength of 2.5 kpc) and a flat outer 
HI velocity dispersion that rises parabolicly in the inner disk 
from an outer plateau of 7 km s"' to a central peak of 10 km s"' . 
In Fig. |4a]we show the XV diagrams for this synthetic galaxy 
without noise and with noise introduced at peak S /N levels of 
80 and 20. The upper diagrams, labelled 'raw unedged', are the 
pure synthetic diagrams that would be observed in the absence 
of noise. The lower ones show the effects of introducing noise. 
The radial (horizontal) scales in the diagrams left and right are 
different and have been chosen such that it the observable radial 
extents fit within the frames. This shows the severe effects of 
noise. 

In Figs. |4b] and |4c] we show the measured kinematics from 
the simulated observation of this model galaxy with noise lev- 
els corresponding to peak signal-to-noise ratios of 80 and 20. In 
each figure the actual kinematics are shown in grey (yellow). As 



mentioned above, both the rotation curve and the radial surface 
density are well measured, even in the low sensitivity S IN- 20 
and S IN- 30 figures. However, the HI velocity dispersion is par- 
ticularly sensitive to high noise levels that submerge the outer 
HI distribution. Comparison of an ideal XV map and XV map 
generated with finite telescope sensitivity and resolution (see 
Fig. |4ali demonstrates the limited HI detected at low sensitivi- 
ties. But despite the low flux recovery, these simulations show 
that low sensitivity observations still yield reasonably accurate 
kinematics, albeit with dramatically increasing measurement er- 
ror at lower sensitivities. 

To reduce the scatter, an evenly weighted running mean was 
applied to the measured fits after the XV diagram was fitted with 
the radial decomposition method. By averaging the measured 
kinematics every 3 or more pixels in radius, it was possible to 
reduce the scatter in measured velocity dispersion and recover 
the radial trend of the kinematics with good confidence. This 
proved to be more revealing, although it effectively lowered the 
spatial resolutions of the measurement. The simulations show 
that measurements of galaxies observed with peak S IN < 100 
benefit from this post-fit averaging. Given an XV image gridded 
with a beam sampling density of 5 (i.e. each pixel is 1 /5 of the 
beam FWHMg) and beamwidth FWHMfl= 300 pc, the final post- 
averaging resolution is still high at 180 or 300 pc after averaging 
over a 3-pixel or 5-pixel wide bin. A post fit averaging interval of 
5 was used for the S IN,„ax- 80 dataset; while this was increased 
to 7 for the noisy S INmax- 20 XV map. As shown by the plots 
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Fig. 8a. Fitted kinematics (with beam correction) measured from XV diagram observed with a telescope beam of FWHMfl= 2 kpc 
and with a peak signal-to-noise of 80. To reduce the scatter in the measured kinematics, the velocity dispersion and rotation curve 
were averaged (post fitting) over an evenly weighted bin of width of 3 pixels. In Fig. [8b] the XV diagram built with the best-fit 
measured kinematics and observed to the same telescope specifications is compared to the original XV diagram. 



in Figs. l4bl and l4cl it was possible to recover the radial trend in 
the velocity dispersion for XV diagrams with these kinematics, 
providing the peak S /N > 20. 

In Figs. l4bl and l4cl we show the measured kinematics from 
the simulated observation of this model galaxy with noise lev- 
els corresponding to peak signal-to-noise ratios of 80 and 20. In 
each figure the actual kinematics are shown in grey (yellow). As 
mentioned above, both the rotation curve and the radial surface 
density are well measured, even in the low sensitivity S IN- 20 
and S IN- 30 figures. However, the HI velocity dispersion is par- 
ticularly sensitive to high noise levels that submerge the outer 
HI distribution. Comparison of an ideal XV map and XV map 
generated with finite telescope sensitivity and resolution (see 
Fig. Hil) demonstrates the limited HI detected at low sensitivi- 
ties. But despite the low flux recovery, these simulations show 
that low sensitivity observations still yield reasonably accurate 
kinematics, albeit with dramatically increasing measurement er- 
ror at lower sensitivities. 

In the slice fits of the edge of the galaxy, the measurements 
are biased by two effects. Firstly, in the outer few slices the radial 



decomposition method models all the flux as originating from a 
single radial bin. Consequently, the measured velocity disper- 
sion is broadened by the projected rotation of gas away from 
the line-of-nodes, i.e. at radii higher than R-R'; similarly the 
measured surface density is overestimated due to the attribution 
of flux away from the lines-of-nodes. This problem is system- 
atic of the radial decomposition method, although it is relatively 
unimportant as it only effects the outer few slices. 

Secondly, the fits to the velocity profile of slices with very 
low peak ( i 5 - 8) are biased by noise peaks dominating 
the underlying flux profile. The radial decomposition fitting al- 
gorithm starts fitting where it finds velocity profile slices with 
flux greater than 2cr in 4 adjacent velocity channels. A galaxy 
with an exponentially declining radial surface density typically 
has a peak S IN i 5 - 8 in the slices at the outermost 1-2 kpc. 
Gaussian fitting of these low S IN velocity profiles is particu- 
larly sensitive to underestimating the velocity dispersion, due to 
fitting the noise peaks instead of the underlying flux profile. As 
noise peaks typically have a FWHM of 1-2 velocity channels, 
the velocity dispersion measurements in these outer few kpc is 
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Fig. 8b. Top: Original XV diagram observed with a FWHMg^ 2 kpc telescope beam and peak S /N of 80. Middle: The XV diagram 
built using the best-fit measured kinematics and observed with the same telescope resolution and sensitivity limits. The measured 
kinematics used to build the model XV plot are shown in Fig. [8a] Bottom: The difference between the top and middle diagrams. 
Contour units are 2, 3, 5, 7, 10, 15, 20, 30, 50, 70 times the rm.s. noise in the channel maps. 



contaminated with noise fits exhibiting a velocity dispersion of 
o",, ~ (1 - 2 channels)/(2 -\/21og2) ~ 1 - 2 km s"'. As the ra- 
dial decomposition fitting method progressively fits inwards us- 
ing the outer fits to model the gas at radii away from the line-of- 
nodes in each velocity profile, erroneous fits cause subsequent 
measurements to oscillate around the correct value. These er- 
roneous fits only effect the adjacent few slice fits as the typi- 
cally near-exponential radial surface density in the outer HI disk 
causes the flux contribution of outer radii to be relatively small. 



In summary, the radial decomposition method is able to suc- 
cessfully measure the radial distributions of the velocity disper- 
sion, rotation velocity and surface density for XV diagram ob- 
servations with peak S /N greater than 20. This could perhaps 
be undertaken with lower S /N maps, if the observation is actu- 
ally detecting the bulk of the actual HI distribution, as could be 
the case on a galaxy with an intrinsically shallow radial surface 
density. The method breaks down for lower S /N maps which are 
dominated by the noise over a larger major axis extent of the de- 
tected gas distribution. Also, as can be seen from the above fig- 
ures, the major limitation of the low signal-to-noise observation 
of the above galaxy model is that one is only detecting the inner 
half of the gas disk, because the outer disk HI lies below the 2cr 
noise level. The outermost radius with sufficient flux above 2cr 
to allow kinematic fitting is marked by vertical dashed lines on 
the XV maps observed with different noise limits in Fig.l4al 



4.2. Effect of spectral resolution 

Modern radio telescope array correlators of the ATCA and VLA 
have a spectral resolution of 3 - 5 km s"', when configured for 
a frequency bandwidth sufficient to cover the rotational velocity 
range of a small-to-medium spiral galaxy (v^ax ~ 200 km s"'). 
This is sufficient to resolve the FWHM due to the typical intrin- 
sic gas velocity dispersion of 7 km s"' (the commonly adopted 
value) by 3 - 5 channels. For this reason we did not investi- 
gate how spectral resolution changes affected the measurement 
accuracy, although it should probably be considered if low HI 
velocity dispersions are found. All measurements of simulated 
HI observations were corrected for instrument broadening, ac- 
cording to Eqn. (|9]l. 

4.3. Effect of spatial resolution 

The shape of the beam is the largest factor apart from the sen- 
sitivity affecting the accuracy of the measured kinematics. But 
the effect of the beam smearing on the measured HI kinemat- 
ics is only serious if the intrinsic galaxy kinematics has a steep 
gradient over the scale of the beam. As a result, the observed 
HI distribution is more distorted over projected radial domains 
with steep gradients of rotation, velocity dispersion or surface 
density than those parts of the galaxy with shallower gradients. 
In practice, the observed HI distribution of a small galaxy with 
a maximum rotation speed of i 100 km s"' and a slowly ris- 
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Fig. 9a. Fitted kinematics (with beam correction) measured from XV diagram observed with a telescope beam of FWHMfl= 2 kpc 
and with a peak signal-to-noise of 30. To reduce the scatter in the measured kinematics, the velocity dispersion and rotation curve 
were averaged (post fitting) over an evenly weighted bin of width of 9 pixels. In Fig. |9b] the XV diagram built with the best-fit 
measured kinematics and observed to the same telescope specifications is compared to the original XV diagram. The plots show 
that the larger beam of this observation allows the shape of the velocity dispersion to be recovered more accurately than a higher 
resolution observation with a 1.0 kpc telescope beam, observed to the same noise limit (compare to Fig. [Tab. 



ing inner velocity dispersion is only minimally distorted by a 
telescope beam with FWHMe= 300 pc. (See Figs. l5alandl5blto 
view the magnitude of the beam-smearing effects of different the 
beam size on the measured kinematics of the model galaxy.) 

Calculation of the dark halo mass density requires smooth 
functions of the kinematics. Consequently, it is more important 
to determine the gas kinematics on scales of 1 kpc, than 100 
pc. The sensitivity of the observations and accuracy of the radial 
kinematics was improved by convolving the observations by a 
larger beam size prior to fitting. This is demonstrated in Figs. |7a| 
and|9a]which show the measured kinematics using beam correc- 
tion from the same galaxy observed with FWHM beams of 1 and 
2 kpc, respectively, at a low sensitivity with peak S /N of 30. 

However, as the beam-smearing correction is performed by 
deconvolution — which amplifies the effective noise — it is im- 
portant to perform beam correction only over the major axis 
range affected by beam smearing. Otherwise one is adding un- 
necessary scatter to the measurements across the whole radial 



extent of the galaxy. To investigate the effect of beam-smearing 
on the measurements kinematics, simulations we analysed sim- 
ulated, noise-less XV maps. Figs.l5alandl5blshow the results of 
fits when beam-smeared with a telescope beam with an FWHMe 
of 1.0 and 2.0 kpc beams. The beam-smearing effect becomes 
quite dramatic in the measurements of galaxies observed with 
large beams (2 kpc or more). Because of the location of the grid 
points, the uncorrected kinematics shows some marked asym- 
metries (e.g. Fig.|5ali which are almost completely corrected by 
the beam-smearing correction process (Figs.l6alandl7ab. 

As explained in Sect. 13.21 beam-smearing causes the veloc- 
ity profile of a slice to deviate from a sum of Gaussians. As a 
result, the measured fits of the velocity dispersion and rotation 
velocity oscillate around the true value, due to the progressive 
nature of the fitting process. Similar oscillations also occur in 
the fitted surface density; however these oscillations are on flux 
scales lower than the scatter due to noise. 
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Fig. 9b. Top: Original XV diagram observed with a FWHMg^ 2 kpc telescope beam and peak S /N of 30. Bottom: The XV diagram 
built using the best-fit measured kinematics and observed with the same telescope resolution and sensitivity limits. The measured 
kinematics used to build the bottom plot are shown in Fig.|9a] Bottom: The difference between the top and middle diagrams. Contour 
units are 2, 3, 5, 7, 10, 15, 20, 30 times the rm.s. noise in the channel maps. 



The best procedure is to first perform an initial radial de- 
composition fit of the galaxy without beam-smearing correction. 
This gives accurate kinematic fits at radii where there is no large 
intrinsic gradients of the kinematics. The radii affected by beam- 
smearing are determined from the radii with oscillating mea- 
sured fits. Subsequent iterations of the fitting process, are per- 
formed using beam smearing in this radial range. As explained 
in Sect. 13.21 the shape of the inner galaxy kinematics are esti- 
mated by extrapolating inwards to build the mini-XV diagram 
used in beam-correcting the radii affected by beam-smearing. 
Fortuitously, the large radial gradients of the intrinsic kinematics 
that necessitate beam smearing correction occur predominantly 
in the inner disk. Most late-type, gas-rich galaxies contain suffi- 
cient HI in the inner disk to apply the beam smearing correction 
method. 

Figs. |6a| to |9a] show the measured kinematics fits obtained 
with beam-smearing correction of XV maps obtained from mod- 
erate (FWHMe= 1.0 kpc) and low spatial resolution (FWHMe= 
2.0 kpc) XV maps. The measurements are undertaken on both 
low and high sensitivity XV maps (peak S IN- 30 and S IN- 80, 
respectively) at each spatial resolution. The XV maps built with 
the measured kinematics are displayed in Figs.l6blto|9bl 

The residual effects of beam-smearing after correction in 
high S/N data (peak S/N - 80) can be appreciated by comparing 
Figs. l4bl(300pc beam). l6al(l kpc beam) and [8a] (2 kpc beam). 
It can be seen that the beam correction process does a good job 
of correcting the beam-smearing that occurs over the inner disk. 



The velocity dispersion is well recovered over large radial scales, 
except for the measurements on the low signal-to-noise {S IN- 
30) dataset observed with the 1 kpc beam. Very high signal-to- 
noise observations (S IN i 100) were able to recover the radial 
velocity dispersion structure on sub-beam scales to accuracies of 
less than 1 km s"', however measurements to observations with 
large beams FWHMe > 1 kpc at 30 i i 100 are limited to 
recovered the velocity dispersion on radial scales of ~1 FWHMg 
and an uncertainty of approximately 1-2 km s"'. 

The measured rotation curve is accurately recovered. The 
surface density is also recovered with good accuracy, though the 
uncertainty is more sensitive to the noise level of the XV map. 

4.4. Fine tuning the iterative method 

To minimise the noise sensitivity of the velocity dispersion fit- 
ting method, various adaptions were trialed on simulated obser- 
vations with different spatial resolutions and noise levels. These 
adaptions trialed were: 

• Binning of the R' dimension of the XV map prior to fitting 
to increase the signal-to-noise. 

• Smoothing of each velocity profile slice through the XV 
map by a Gaussian of known dispersion prior to fitting kinemat- 
ics to each velocity slice. Profile broadening by a Gaussian with 
varying FWHM from 3 to 6 velocity channels was tried. This 
additional known source of velocity dispersion was taken into 
account in the derivation of the HI velocity dispersion. 
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Fig. 10a. These plots display measures of the fit to the kinematics and how these are affected by varying the iterative onion peeler 
for a synthetic galaxy with signal-to-noise S /N= 80, and a small beam Osy„= 0.3 kpc. The first three columns show the integrated 
absolute difference between the final model and the the original XV diagram for each kinematic function (from left to right respec- 
tively velocity dispersion 'SigmaR', rotation velocity ' V and surface density 'Surf'), normalised by the number of sightlines within 
a fit. The units for the total fit difference of the velocity dispersion and the rotation curve are in km s ', while the surface density 
difference is measured in log flux. The fourth column displays the integrated flux difference between the final XV diagram built 
from the best fit kinematics, and the original XV diagram, as a percentage of that in the original XV diagram. See the text for further 
explanation. 



• Smoothing of the derived fit cr^.HiiR) with a Gaussian 
FWHM of varying number of R' pixels. The smoothed fits were 
then compared to determine the new model o-y HiiR)- Only the 
velocity dispersion fit was smoothed, as the rotation curve and 
deprojected surface density fits exhibited a very low scatter. 



• Damping of the model incrementations A (see Fig. O to 
minimise oscillations in the fits used to create the new model. 
Damping was only applied to the model incrementation of the 
velocity dispersion, as the small scatter of the fits to the rotation 
curve and HI surface density removed the need for any model 
damping. 
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Fig. 10b. These plots display the same measures of the fit to the kinematics as in Fig. llOal but now for a synthetic galaxy with 
signal-to-noise S /N -30. Lay-out and units are the same as in Fig. llOal 



Several variations were tried for each adaption, resulting in 
108 separate simulations, which were run on both a high (S IN- 
80) and low (S IN- 30) signal-to-noise XV map observed with 
a 0.3 kpc beam. The galaxy model used was a exponential HI 
surface density, constant HI velocity dispersion and a differential 
rotation curve. Two indicators of the accuracy were measured. 

Firstly, a residual XV diagram was calculated by subtracting 
the original XV diagram from the XV diagram formed from the 
best fit kinematics. The summed absolute residual XV flux in % 
of the total flux in the maps is displayed in the plots in the far 
right column of Fig. fTOal for the S IN- 80 results, and Fig. llObl for 
the S IN- 30 results. While the residual XV diagram is a good 



indicator of measurement accuracy on large scales, it is quite in- 
sensitive to measurement errors on small radial scales. Except 
for one fit which diverged in the simulations on low signal-to- 
noise XV diagram, the summed residual was not more than 1 .0% 
(0.25% in the SN80 dataset) of the total flux in the original map 
for both the low and high S IN measurements, being approxi- 
mately three times greater in the fit to the low signal-to-noise 
map. 

Secondly, we also used the difference between the kinemat- 
ics measured from the XV diagram formed with the best fit kine- 
matics and the kinematics measured from the observed XV dia- 
gram. The difference in fit (again absolute values) at each radius 
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was then summed to form a total difference for the fit to that 
kinematic function. We were able to compare the fits obtained 
from binned and un-binned XV maps, by normalising the total 
difference in fit by the number of sightlines used in the fit. The 
normalised total fit difference for each kinematic function is dis- 
played in the plots in the first three columns of Fig. llOal for the 
S /N= 80 results, and Fig. [TOb]for the S /N= 30 results. 

Each line marks the trend in the normalised total fit differ- 
ence if only the method plotted on the abscissa is varied, and 
all other method variations are held constant. Sometimes the 
diamond symbols overlap, giving the impression that there are 
fewer simulations in some of the panels. However, the number 
of lines in each panel is the same and the variation of parameters 
is distributed evenly over the various panels. The large spread 
in normalised total fit difference of the different lines is primar- 
ily due to binning the XV diagram along the major axis prior 
to fitting. In the plots showing how the measurements varied 
with other method variations, the lines with high measurement 
error are those which used a large major axis binning interval. 
As expected, binning the XV diagram prior to fitting reduced 
the amount of major axis structure information, which dramat- 
ically increased the measurement uncertainty of all three fitted 
functions. The different values of the binning interval have in 
some cases resulted in the outcomes being distributeded in two 
somewhat distinct groups. 

The plots show that smoothing the fits before incrementing, 
and damping the model incrementation offer no improvement. 
For the high signal-to-noise dataset, the velocity dispersion mea- 
surements were improved by smoothing the velocity profile of 
each slice by a small Gaussian prior to fitting. The velocity dis- 
persion measurements to the peak S /N- 80 dataset were im- 
proved by 25% (or on average 0.5 - 1 .0 km s"' ) with a Gaussian 
smoothing kernel with a FWHM of 2 velocity channels, i.e. 
slightly larger than the FWHM of noise which will typically be 
1-2 channels. However, Gaussian smoothing the velocity axis 
caused no distinct improvement in the measurements to the 
low signal-to-noise dataset. This is probably because Gaussian 
smoothing a noisy profile distorts the shape slightly by incorpo- 
rating the noise into the smoothed profile shape. Consequently, 
we have not adopted any additional smoothing or other varia- 
tions discussed in this section (Sect. 14.41 ) to the radial decompo- 
sition method. 



4.5. Accuracy of the iterative radial decomposition metinod 

Convergence to the best kinematic model was usually achieved 
after three iterations of the radial decomposition fitting process. 
The average measurement errors of each kinematic fit deter- 
mined with the adopted unchanged radial decomposition method 
are: Act,. h/= 0.7, 1 .3 km s^' (for peak S /N^ 80 and peak S /N^ 
30, respectively); AV^ 0.9, 1.4 km s"', and A(log !///)= 0.1, 
0.15. The sum of the residual XV diagram shows that the ra- 
dial decomposition level was able to recover the flux to a level 
of 0.2% and 0.6% of the total flux of the observed XV diagram 
(for maps with peak S /N= 80 and S /N- 30, respectively). This 
measurement accuracy is indicative only, as it constitutes a mean 
error to measurements of a galaxy with this particular kinematics 
and structure. As mentioned earlier, fits to a XV map observed 
to the same sensitivity level of a galaxy with a less steep intrinsic 
surface density should have a lower measurement error. 



5. Summary 

Using our new iterative radial decomposition method, we show 
that it is possible to accurately measure all three functions that 
describe the planar kinematics and distribution of a galaxy: the 
radial HI surface density, the rotation curve and the HI velocity 
dispersion. This method is a significant improvement on previ- 
ous methods to model edge-on galaxies which approximated one 
or more of these kinematic functions. Measurements of simu- 
lated galaxies show that a peak signal-to-noise of i 30 is needed 
to accurately measure all three kinematic functions on a galaxy 
with a intrinsically steep (exponential over all R) HI surface den- 
sity. However, reliable measurements would probably also be 
obtained with observations at a peak signal-to-noise (S /N i 20) 
Furthermore, the method should also be successful on galaxies 
with an intrinsically shallower surface density — as is common 
in most disk galaxies due to depletion of HI in the stellar disk. 

In the next paper in this series, we present the derived HI 
kinematics and radial surface density of our galaxies as mea- 
sured with this method. Paper III also contains the vertical HI 
flaring results, which is measured as a function of galactocentric 
radius using the deprojected HI distribution and kinematics. 
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